Nuclear Superfluidity and Specific Heat in the Inner Crust of Neutron Stars 
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We analyse the temperature dependence of pairing correlations in the inner crust matter of neutron 
stars. The study is done in a finite-temperature HFB approach and by using a zero range pairing 
force adjusted to the pairing properties of infinite neutron matter. Within the same approach we 
investigate how the specific heat of the inner crust depends on temperature, matter inhomogeneity, 
and the assumption used for the pairing force. It is shown that in a physical relevant range of 
densities the pairing properties of inner crust matter depend significantly on temperature. The finite- 
temperature HFB calculations show also that the specific heat is rather sensitive to the presence of 
nuclear clusters inside the inner crust. However, the most dramatic change of the specific heat is 
determined by the scenario used for the neutron matter superfluidity. 



The inner crust of neutron stars consists of a lattice of 
neutron-rich nuclei immersed in a sea of unbound neu- 
trons and relativistic electrons. Down to the inner edge 
of the crust, the crystal lattice is most probably formed 
by spherical nuclei. More inside the star, before the nu- 
clei dissolve completely into the liquid of the core, the 
nuclear matter can develop other exotic configurations 
as well, i.e, rods, plates, tubes, and bubbles (for a discus- 
sion of inner crust matter see p] and references therein). 
The thickness of the inner crust is rather small, of the 
order of one kilometer, and its mass is only about 1% 
of the neutron star mass. However, in spite of its small 
size, the properties of inner crust matter, especially its 
superfluid properties, have important consequences for 
the dynamics and the thermodynamics of neutron stars. 

The superfluid properties of the inner crust have been 
considered long ago in connection withpost-glitch timing 
observations and colling processes @,y,|j] However, al- 
though the neutron star matter superfluidity has been in- 
tensively studied in the last decades [f| , so far only a few 
microscopic calculations have been done for the superflu- 
idity of inner crust matter. The existing calculations are 
done in the framework of the Hartree-Fock-Bogoliubov 
(HFB) approach and at zero temperature |E0H)I3- As 
it is well-known, the limitation to zero temperature is 
justified only if the pairing gap of the superfluid is much 
larger than its temperature. This is not generally the 
case for the inner crust matter, where this condition is 
fulfiled only for some limited density regions. 

The purpose of this paper is to study how the pair- 
ing properties of inner crust matter are changing with 
the temperature and what are the consequences of these 
changes on the specific heat. The present investigation 
is done in a finite-temperature HFB approach. In the 
calculations both the mean field and the pairing field are 
determined self-consistently by using a Skyrme-type in- 
teraction and an effective pairing force adjusted to the 
superfluid properties of infinite neutron matter. The set- 
ting of the calculations follows Ref.0, where the same 
two-body forces have been used for describing the pair- 
ing properties of inner crust matter at zero temperature. 

The finite-temperature HFB (FT-HFB) equations have 



a structure similar to the HFB equations at zero temper- 
ature ,10]. For zero range forces and spherical symmetry 
the radial FT-HFB equations can be written as: 
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where Ej is the quasiparticle energy, Ui, Vf are the com- 
ponents of the radial FT-HFB wave function and A is 
the chemical potential. The quantity hx(r) is the ther- 
mal averaged mean field hamiltonian and Ay(r) is the 
thermal averaged pairing field, fn a self-consistent cal- 
culation based on a Skyrme-type force as used in this 
study, hr(r) is expressed in terms of thermal averaged 
densities, i.e., kinetic energy density r-r(r), particle den- 
sity pt{t) and spin density Jr(r), in the same way as 
in the Skyrme-HF approach If]. The thermal averaged 
densities are given by: 
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where fi = [1 + exp(Ei / ksT)] -1 is the Fermi distribu- 
tion, kg is the Boltzmann constant and T is the tem- 
perature. The summations in the equations above are 
over the whole positive energy quasiparticle spectrum, 
including the continuum states. In the latter case the 
summations should be replaced by an integral over the 
energy 
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The thermal averaged pairing field is calculated with 
a density dependent contact force of the following form 

m 

V(r-r>) = Voil-vi^-mr-r') = V eff (p(r))S(r-r>), 
Po 

where p(r) is the baryonic density. With this force the 
thermal averaged pairing held is local and is given by: 

A T (r) = V eff (j>(r))±'£(2j i + l)Ut(r)V i (r)0.-2f i ) 

i 

ee V eff { P {r))K T {r), (6) 

where kt(^) is the thermal averaged pairing tensor. Due 
to the density dependence of the pairing force, the ther- 
mal averaged mean field hamiltonian hx{r) depends also 
on kt- 

The inner crust matter is usually divided into indepen- 
dent Wigner-Seitz cells, each cell containing in its cen- 
ter a neutron-rich nucleus surrounded by unbound neu- 
trons and permeated by a relativistic electron gas uni- 
formly distributed inside the cell. In this study the FT- 
HFB equations are solved for some representative spher- 
ical cells determined in Ref . [T3| . To generate far from 
the nucleus a constant density corresponding to the neu- 
tron gas, the FT-HFB equations are solved by imposing 
Dirichlet-Neumann boundary conditions at the edge of 
the cell 0] , i.e., all wave functions of even parity vanish 
and the derivatives of odd-parity wave functions vanish. 

In the FT-HFB calculations we use for the particle- 
hole channel the Skyrme effective interaction SLy4 [Tlj . 
which has been adjusted to describe properly the mean 
field properties of neutron-rich nuclei and infinite neu- 
tron matter. For the effective pairing force the choice 
is more problematic. A realistic FT-HFB calculation of 
crust superfluidity should be based on a pairing force able 
to describe reasonably well at least the pairing properties 
of neutron matter, which forms the main part of the inner 
crust. However, the magnitude of pairing correlations in 
neutron matter is still a subject of debate. On one hand, 
BCS calculations based on a Gogny force commonly used 
in finite nuclei ^(|, give for the pairing gap of infinite 
neutron matter a maximum value of about 3.2 MeV at 
a Fermi momentum kp w 0.9 fm -1 17]. On the other 
hand, the microscopic calculations based on induced in- 
teractions predict for the maximum value of the gap much 
smaller values, of around 1 MeV 0,0,01- These two 
different scenarios for the neutron matter superfluidity 
are used here to determine the parameters of the pairing 
force (5) employed in the FT-HFB calculations. Thus, for 
each Wigner-Seitz cell we perform two FT-HFB calcula- 
tions, with two sets of parameters for the pairing force. 
For the first calculation we use the parameters: Vo=-430 
MeV fm 3 , rj—0.7, and a=0.45. With these parameters 
and with a cut-off energy for the quasiparticle spectrum 
equal to 60 MeV one obtains approximately the pa i ring 
gap given by the Gogny force in nuclear matter |13l |20| . 



In the second calculation we reduce the strength of the 
force to the value Vo=-330 MeV fm 3 . With this value of 
the strength we simulate the second scenario for the neu- 
tron matter superfluidity, in which the maximum gap in 
neutron matter is around 1 MeV. In principle, the other 
two parameters of the pairing force should be adjusted 
to the density dependence of the gap in neutron matter 
corresponding to the calculations based on induced in- 
teractions. However, since at present this dependence is 
not well-established, we keep for the parameters rj and a 
the same values as in the first calculation. 

The FT-HFB results are analysed here for two repre- 
sentative Wigner-Seitz cells chosen from Ref.0]. These 
cells contain Z=50 protons and have rather different 
baryonic densities, i.e., 0.0204 fm~ 3 and 0.00373 fm~ 3 . 
The cells, which contain N=1750 and N=900 neutrons, 
respectively, are denoted below as a nucleus with Z pro- 
tons and N neutrons, i.e., 1800 Sn and 950 Sn. For the bary- 
onic densities and the nucleonic numbers listed above, 
the radii of the cells are equal to 27.6 fm and 39.3 fm, re- 
spectively. In the calculations we integrate the FT-HFB 
equations up to 27.6 fm for both cells. Up to this dis- 
tance in the second cell there are N=279 neutrons. Since 
at 27.6 fm we are already far out in the region of the 
neutron gas, the pairing field of the second cell is not 
influenced by this choice of the integration radius. 

The FT-HFB calculations are done up to a maximum 
temperature of T=0.5 MeV, which is covering the tem- 
perature range of physical interest. According to numeri- 
cal simulations |2l| , the temperature in the inner crust is 
increasing from about T=0.1 MeV in the low-density re- 
gion to about T=0.3 MeV at the inner edge of the crust. 
This temperature profile corresponds to the state of mat- 
ter before the cooling wave is passing through the inner 
crust from the core region to the star surface. 

The temperature dependence of the pairing fields in 
the two cells presented above is shown in Figs. 1-2. First, 
we notice that for all temperatures the nuclear clusters 
modify significantly the profile of the pairing field. Thus, 
for the high density cell 1800 Sn we can see that the pair- 
ing field is strongly decreasing (in absolute value) in the 
cluster region. However, as already noticed in zero- tem- 
perature HFB calculations , the situation is completely 
different for the low-density cell 950 Sn, where the pair- 
ing field is strongly increasing in the transition region 
from the neutron gas to the nuclear cluster. These differ- 
ent behaviours of the pairing fields are connected to the 
density dependence of the pairing gap in neutron matter 

0E1GI. 

As expected, the dependence on temperature of the 
pairing field is weak when the gap values are much larger 
than the temperature. This is the case for the high- 
density cell 1800 Sn and the first pairing force. For all 
other cases the temperature dependence of the pairing 
field is significant. This is clearly seen in the low-density 
cell 950 Sn. From Fig.l we can see that for the second pair- 
ing force even in the high-density cell 1800 Sn the pairing 
field is changing significantly with the temperature. 
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FIG. 1: Neutron pairing fields for the cell 1800 Sn calculated 
at various temperatures. The numbers 1 and 2 which fol- 
low the cell symbol (see the inset) indicate the variant of the 
pairing force used in the calculations. The full and the dashed 
lines corresponds (from bottom upwards) to the set of tem- 
peratures T={0.0, 0.5} MeV and T={0.0, 0.1, 0.3, 0.5}MeV, 
respectively. 
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FIG. 2: The same as in Fig.l, but for the cell 950 Sn. The 
full and the dashed lines corresponds (from bottom up- 
wards) to the set of temperatures T={0.0, 0.1, 0.3, 0.5} and 
T={0.0, 0.1}MeV, respectively 



To analyse how the pairing correlations are globally 
changing in the cell as a function of temperature, in Fig. 3 
we plot the pairing energy per neutron for the cases in 
which its temperature dependence is significant. In Fig. 3 
are shown also the results obtained when in the same 
cell are distributed uniformly only the neutrons. One 
can notice that the pairing energy of the non-uniform 
system, i.e., cluster plus the neutron gas, has a rather dif- 
ferent temperature dependence compared to the uniform 
neutron matter. From Fig. 3 we can eventually extract 
a critical temperature for the superfluid-normal phase 
transition in the non-uniform nuclear matter of the cell. 



However, as seen for the cell 950 Sn, we should be aware of 
the fact that close to this averaged critical temperature 
the non-uniform system may still have regions in which 
the pairing field has non-negligible values. 
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FIG. 3: Pairing energy per neutron as a function of temper- 
ature. The calculated values, done in step of AT=0.1 MeV, 
are joined by straight lines. The symbol "n" following the 
cell symbol (see the inset) refers to the results obtained when 
in the corresponding cell are distributed uniformly only the 
neutrons. The rest of the notations are the same as in Fig.l. 



We turn now to the specific heat provided by the FT- 
HFB calculations. The specific heat of a given cell of 
volume V is defined by: 



C v = 



1 d£(T) 
V dT ' 



(7) 



where £(T) is the total energy of the baryonic matter 
inside the cell, i.e., 



£(?) = J2 fi E i 



(8) 



Due to the energy gap in the excitation spectrum, the 
specific heat of a superfluid system is dramatically re- 
duced compared to its value in the normal phase. Since 
the specific heat depends exponentially on the energy 
gap, its value for a Wigner-Seitz cell is very sensitive 
to the local variations of the pairing field induced by the 
nuclear clusters. This can be clearly seen in Fig. 4, where 
the specific heat is plotted for the cell 1800 Sn and for the 
neutrons uniformly distributed in the same cell. One can 
notice that at T=0.1 MeV and for the first pairing force 
the presence of the cluster increases the specific heat by 
about 6 times compared to the value for the uniform neu- 
tron gas. A comparable effect of the nuclear cluster on 
the specific heat of the cell 1800 Sn was found earlier by 
using in Eq. (8) the quasiparticle spectrum provided by a 
zero-temperature HFB calculation 0, Q ■ This approxi- 
mation is justified only when the temperature variation 
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of the pairing field in the Wigner-Seitz cells is small, as 
in the case of high density cells and a Gogny type force. 
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FIG. 4: Specific heat for the cell 800 Sn as a function of 
temperature. The notations used in the inset and the repre- 
sentation of the calculated values are the same as in Figs. 1-3. 



The behaviour of the specific heat for the low- density 
cell 950 Sn is shown in Fig. 5. For the first pairing force 
we can also see that at T=0.1 MeV the cluster increases 
the specific heat by about the same factor as in the cell 
i800g n However, for the second pairing force the situ- 
ation is opposite: the presence of the nucleus decreases 
the specific heat instead of increasing it. From Fig. 5 we 
can also see that the specific heat predicted by the two 
pairing forces are not very much different (notice the ex- 
panded scale) for this low-density cell. Due to the fast 
decreasing of the pairing correlations with the tempera- 
ture, at T=0.3 MeV all calculations provide almost the 
same specific heat. 

We have also performed FT-HFB calculations for two 
other cells with 40 protons and having the densities equal 
to 0.0475 fm- 3 and 0.00159 fm" 3 . We found that the 
temperature dependence of the pairing field and specific 
heat in the first and in the second cell with Z=40 is very 
similar to the one in the cells 1800 Sn and 950 Sn, respec- 
tively. 
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FIG. 5: The same as in Fig.4, but for the cell 950 Sn. 

The most striking fact seen in Fig.4 is the huge differ- 
ence between the predictions of the two pairing forces. 
Thus, for T=0.1 MeV this difference amounts to about 7 
orders of magnitude. This discrepancy between the pre- 
dicted values of specific heat, generated by the changes 
in the pairing gaps, is decreasing rapidly with the tem- 
perature, but it remains significant up to T=0.5 MeV. 
From Fig.4 we can see also that for the cell 1800 Sn the 
effect of the cluster on specific heat is similar for both 
pairing forces. 



In conclusion, we have shown that the superfluid prop- 
erties and the specific heat of inner crust matter are af- 
fected rather strongly and in a non-trivial way by the 
temperature of the crust. Therefore, the temperature de- 
pendence should be included in the microscopic calcula- 
tions aiming at describing accurately the superfluid prop- 
erties of inner crust matter. However, the most serious 
problem in microscopic calculations remains the pairing 
interaction. As seen above, shifting to the second pair- 
ing force adjusted to the nuclear matter properties calcu- 
lated with induced interactions, the specific heat of high 
density cells increases by several orders of magnitude. 
Consequently, for the second pairing force the neutron 
and the electron contributions to the total specific heat 
of high-density cells are becoming competitive, at vari- 
ance with the results based on a Gogny force scenario 
8J. To analyse the implications of these two different 
scenarios for the pairing force on physical properties as 
the cooling time one needs to calculate within the FT- 
HFB approach the heat diffusion along the whole inner 
crust. The succesive step would be to include also into 
the microscopic calculations the ground state correlations 
of the inner crust matter, which can modify significantly 
the neutron single-particle properties related to the heat 
diffusion (e.g., effective mass and the mean free path) 
[2^. It is to be expected that such microscopic calcula- 
tions could be eventually used to determine accurately 
the consequences of inner crust superfluidity on the cool- 
ing time of neutron stars. 
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